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Abstract — Rate outage probability is an important performance 
metric to measure the level of quality of service (QoS) in the 
Generation (4G) broadband access networks. Thus, in this 
paper, we calculate a closed form expression of the rate outage 
probability for a given user in a down-link multi-hop OFDMA- 
based system encountered as a result of links' channel 
variations. The channel random behavior on different 
subcarriers allocated to a given user is assumed to follow 
independent non-identical Nakagami-m distributions. Besides 
the rate outage probability formulas for single hop and multi- 
hop networks, we also derive a novel closed form formulas for 
the moment generating function, probability distribution 
function (pdf), and the cumulative distribution function (cdf) of 
a product of independent non-identical Gamma distributed 
random variables (RVs). These RVs are functions of the 
attainable signal-to-noise power ratio (SNR) on the allocated 
group of subcarriers. For single-hop scenario, inspired by the 
rate outage probability closed formula, we formulate an 
optimization problem in which we allocate subcarriers to users 
such that the total transmission rate is maximized while 
catering for fairness for all users. In the proposed formulation, 
fairness is considered by guaranteeing a minimum rate outage 
probability for each admitted user. (Abstract) 

Keywords-component; OFDMA; Nakagami-m; Outage 
probability; Moment generating function 

I. Introduction 

The Worldwide Interoperability for Microwave Access 
(WiMAX) and Long-Term Evolution (LTE)- also known as 
Evolved Universal Terrestrial Radio Access (E-UTRA)- are 
two promising candidates for providing ubiquitous IP 4G 
wireless broadband access. The LTE and WiMAX are 
designed as all IP packet switched networks in order to be 
capable of providing higher link data bit rate [1], [2]. Thus, 
all services involving real and non- real time data, are 
implemented on top of the IP as a network layer protocol. 
Real time data services transmission such as voice and video 
require maintaining the wireless link data rate within the 
requirement of the rate demands of these applications. Thus, 
the outage of the link data rate more than a predefined 
threshold is not tolerable. The link data rate is said to be in 
an outage if the rate requirements of an application cannot 
be satisfied. Hence, we define the rate outage probability as 
the probability of having the attainable link data rate below 



the required application data rate. LTE and WiMAX both 
implement OFDMA as the technique to share the access of 
the media among users. 

OFDMA is a multiple-access scheme that inherits the 
ability of Orthogonal Frequency Division Multiplexing 
(OFDM) to combat the inter- symbol interference and 
frequency selective fading channels. In order to reduce the 
effect of selective fading channels in wireless networks, 
OFDMA divides the total frequency bandwidth into a 
number of narrow-band subcarriers in such a way that each 
subcarrier (SC) experiences a flat fading channel. This is 
possible by selecting the SC width to be close to the 
coherence bandwidth of the channel. A small end-portion of 
the symbol time, which has duration equal to the maximum 
delay spread channel, is called a cyclic prefix (CP). The CP 
is copied at the beginning of the symbol time duration to 
eliminate the effect of inter- symbol interference (ISI). Thus, 
the inter- symbol interference is dealt with while 
orthogonality among the subcarriers is maintained. 

Accordingly, OFDMA is a preeminent access technology 
for broadband wireless data networks compared with 
traditional access technologies such as TDMA and CDMA. 
The main advantages of OFDMA over TDMA/CDMA stem 
from the scalability of OFDMA, the uplink orthogonality 
and the ability to take advantage of the frequency selectivity 
of the channel. Other advantages of OFDMA include 
MIMO-friendliness and ability to support prominent quality 
of service (QoS). 

The rate outage probability defined as the probability of 
having the attainable link data rate below the required 
application data rate. Our objective in this work is to derive a 
closed form expression of the rate outage probability for 
Orthogonal Frequency Division Multiple Access (OFDMA) - 
based system under the Nakagami-m channel model. Our 
interest in calculating the rate outage probability in an 
OFDMA system is motivated by the fact that OFDMA is the 
common physical layer technology for the 4th generation 
(4G) and beyond broadband wireless networks; mainly LTE- 
Advanced and IEEE802.16m networks. Calculating a closed- 
form rate outage probability expression facilitates designing 
sound radio resource management schemes for the efficient 
operation of these networks in single and multi-hop contexts. 

Recently, several attempts are made to compute the 
outage probability in OFDMA systems. These attempts can 
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be categorized into two major categories: the simulation 
modeling and the analytical derivation. An example on 
simulation modeling is the work presented in [3]. Simulation 
requires collecting sufficient statistics of the network, which 
is time consuming, then applying a prediction model, which 
risks various prediction error. Hence, we will focus on the 
analytical derivation category. The works in [4, 5, 6, 7, 8, 
9, and 10] are examples of the analytical derivation attempts 
for single and multi-hop networks. The authors of [4] and [5] 



compute uplink and downlink outage probabilities in [4] and 
[5] respectively under Rayleigh fading channel model. The 
proposal presented in [6] computes downlink outage 
probability of all users in a cell for a particular resource 
allocation (subcarrier). In [7], the author computes the 
downlink outage probability for a system with a limited 
channel feedback, and the proposal presented in [8] attempt 
to describe downlink outage probability related to bandwidth 
user demand. The work in [10], and [9] addressed the 
calculation of outage probability in multi-hop networks. The 
authors of [10] provided three lower bounds of end-to-end 
outage capacity over Rayleigh fading channels, while [9] 
provided an approximation on outage probability through the 
numerical inversion of a Laplace transform over Rayleigh 
fading channels as well . 

The aforementioned analytical proposals, although 
compute outage probability for different objectives, they do 
share common procedures. They compute the outage 
probability for a particular resource allocation (OFDMA 
subcarrier) and/or employed a simplified channel model 
such as the Rayleigh channel model. Moreover, all proposals 
are semi-analytical where outage probability computation is 
approximated through calculating upper and/or lower 
bounds or through providing numerical solutions. 
Nevertheless, none of the above proposals provides a closed 
form expression of the outage probability. In this paper, we 
address the void of the current literature by deriving a closed 
form formula for the rate outage probability in single (i.e., 
PMP) and multi-hop networks. We compute the outage 
probability based on a generalized channel model, namely 
the Nakagami-m model of which Rayleigh is merely a 
special case. Finally, we express the outage probability 
formula as a direct function of the user required data rate and 
two other parameters; the Nakagami-m fading coefficient 
and the spreading factor. 

Our derivation has several applications in 4G and beyond 
broadband networks. For example, outage probability can be 
considered as a QoS metric to meet a specific connection 
data rate requirement, or it can be utilized as a performance 
measure to evaluate the level of meeting the total demands 
of users in a cellular network. 

The remainder of this paper is organized as follows. In 
Section 2 we present the problem description and system 
model, while the computation of the outage probability and 
the relevant mathematical formulation is offered in Sections 
3 and 4. In section 5, we introduce an optimization 



formulation for single hop cellular system scenario. In 
Section 6 we conclude the paper. 

II. SYSTEM MODEL 

We consider an OFDMA system model in two network 

scenarios: namely, single hop cellular radio network, see 

Fig. 1 and multi-hop cellular radio network, see Fig. 2. Each 

cell, the network consists of a Base Station (BS) with 

omnidirectional antenna located in the center of the cell 

serving multiple Subscriber (SS) stations uniformly 

distributed over the cell. 

The SSs share the total channel bandwidth, which is 

divided into n orthogonal subcarrier s. SSs and BS 

transmission is on a time frame basis, where each frame 

consists of multiple OFDMA symbols. The channel 

condition of the SS is assumed not changing during each 

time frame. However, it may change from one frame to 

another. The received low-pass signal at SS k on the nth SC 

at OFDMA symbol t is denoted by ^ ^ and given as: 



%n< r \ 



R k ,n(t) = K n (t)e "*■» S Kn (t) + z Kn (t) 



(l) 



where S kn (t) is the equivalent low-pass transmitted signal 

to the kth SS on the nth SC. h kn (t)e ~^* (0 is the time- 
variant transfer function of the equivalent low-pass 
frequency non- selective channel of the k th SS experienced 

on the nth SC. The transfer function n kn \t)e ' is a 

complex- valued Gaussian random process in which h k n (t) 

is a Nakagami-m distributed random variable and </> kn (t) is 

uniformly distributed random variable in {—71, 7l) . h kn (t) 

represents the channel envelope (modulus), while (/> kn {t) 

captures the phase of the equivalent low-pass channel. We 
Adopt the Nakagami-m distribution because it fits the 
realistic statistical variations of the channel envelope in both 
line-of-sight (LOS) and non line-of-sight (NLOS) 
communications. Both, the envelope and the phase of the 
channel are assumed to be jointly independent. Furthermore, 

we assume that h kn (t) and h km (t) are independent 

random variables Vm ^ n , which implies that the channel 
coherence bandwidth is approximately equal to the 

bandwidth of a SC. Z kn (t) is the complex-valued white 

Gaussian noise process with power spectral density denoted 

by N n • 
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We denote to the SS's k , (ks{l,...,K} on subcarrier n , 
n£{l,....,N}) data rate by r k (t) bits per OFDMA 
symbol. r kn {t) depends on the allocated power p kn , SS 
k on subcarrier n , and the channel gain h k n . The channel 

gain h kn is modeled as a Nakagami-m random variable 
with (a probability density function) pdf defined as: 



m, , m k,n , 2 

1™ k,n ( —)(o 

Zm h ' 2m, -1 Q 

A, (*>) = — ^ <? 



"M 



*.« -'%«-„ "*■« ,VneN,(2) 



rK,„)£V„ 



where Q^ w =£'{/z^ w } is a coefficient that reflects the 
spread of the Nakagami-m distribution and 
N = { 1 , 2, • • • , N} is the indexing set of the data 



subcarriers. m k n is the fading coefficient that measures the 
severity fading of the channel. m k n is defined as: 



Q 



kji 



E{(hl n -n k j 2 } 



(3) 



Note that for m k = oo , this corresponds to a non- 
fading channel. For m k = 1 , the Nakagami-m distribution 

reduces to the Rayleigh distribution, which is used to model 
small scale channel variation in NLOS systems, and to one- 



sided Gaussian distribution for 



. Further 



examples on the generality of the Nakagami-m model is 
when m k n < 1 , in this case Nakagami-m distribution 

closely approximates the Hoyt distribution, and for m k n > 1 

, it approximates the Rician channel model. Thus, the 
Nakagami-m channel model is capable of capturing a 

realistic wireless channel through the two variables; Jfl k n 

and Cl k n . Hence, Nakagami-m channel is an appropriate 

model to meet this work objective, which is to facilitate 
practical resource management through trackable closed- 
form expression of the rate outage probability. 



III. Derivation of Single-Hop 
Rate Outage Probability 

In this paper, the measure of the level of QoS 
satisfaction is given in terms of the attainable data rate , r k (t) 

, against the SS k minimum required data rate, r k min . 

Explicitly, a SS is in outage, if the required data rate is not 
met due to low attainable data rate. The outage of meeting the 
required data rate can be considered as a QoS metric or as a 
performance measure. For the case of the outage probability 
as a QoS metric, a SS can demand a predefined threshold for 
the probability of the outage of service. The outage 
probability threshold will be used to control the amount of 
network resources allocated for the user. Alternatively, the 
outage probability can be calculated given the allocated 
resources. In this case the outage probability is considered as 
a performance measure. For both cases, our aim is to 
calculate the rate outage probability for a specific SS in 
consequence of its variable channel quality and adaptive 
transmission. 

The rate outage probability P™ of the kth SS can 
be expressed as: 

Pr = PKr k (t)<r kMn ) (4) 

Now, let us define the achievable transmission data rate by 
subcarrier n of SS k as [11]: 



r k , n (t) = — \og 2 (l + y kA (t)), 



(5) 



where B is the link bandwidth and N is the total number 
of data subcarriers. Total link bandwidth is denoted as B 

B 

and hence, — represents the physical bandwidth of each 

subcarrier. y k n (t) is the Signal to Noise power Ratio (SNR) 

of SS k at subcarrier n at OFDMA symbol t and is 
expressed as: 

n * {t)= ^JTT> (6) 

where p kn is the transmit power of the SS k at the 

subcarrier n, H k n (t) @ h k n (t) with h k n (t) the channel 

gain between the BS and SS k at subcarrier n at OFDMA 

symbol t . N o is the Additive White Gaussian Noise 

(AWGN) power spectral density at the receiver. As per our 
earlier assumption that the channel gains on all subcarriers 
change form one time frame to another time frame, we drop 
the time index t in our next calculations. Note that in (6), the 
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interference coming from other users using same subcarriers 
as SS k in the same OFDM A symbol has been assumed 
negligible. The impact of the interference in OFDM A 
signalling can be taken care of in the MAC protocol. 

The overall achievable transmission data rate, r k in 

its OFDM A symbol, of SS k on all its allocated subcarriers 

in the subset Q : fl = l 9 ....M k ',M k < N is: 



nsq 



B 



= T7Z lo g2( 1+ n,J> 



N 



(7) 
(8) 



nsC k 



= — log 2 



where M k =\ Q I , i.e. the cardinality of the subset Q . 

Using (9) and (4) the rate outage probability P™ can be 
expressed as: 



pout = pr 



K nGc k 



, (10) 



where B 
Define 



B 



N 



is the bandwidth of each subcarrier. 



(r IB \ 

Y k \c k <2 k > min SC V (15) 



Y k is a product of independent not necessarily identically 

(i.n.i.d) Gamma distributed random variables, X k n . Thus, to 
calculate the rate outage probability, it is required to find the 
probability density function (pdf) of Y k for a given allocated 

subcarriers' subset Q . To derive the pdf of Y k , we first 

derive the moment generating function (MGF) A Y (s) of 

k 

the random variable Y k as illustrated by the following 
theorem: 

Theorem 1: The MGF of Y k defined as a product of i.n.i.d 
Gamma distributed random variables is given by 
1 



Av(j) = 



n r K,j 

n & C k 



(16) 



xG M J 



A,M 



-4IA, 



l-m k l ,...,l-m i 



l k,M, 



where G M \ [.] is a Meijer G-function which is included as 

a standard built-in function in most of the mathematical 
software packages such as MATHEMATICA and MAPLE. 



( U ) Proof: The MGF A Y (s) is the expectation of e k , i.e. 



Given that y k n is based on a Nakagami-m channel model, it 
follows that the random variable X kn has the following 



Gamma distribution: 



m, -1 

jc-r 



fx k W = ^— L exp(— — ), (12) 



where 



1 1 Plr O. 

A = E{X k J = (1 + %^-). 



m ; 



kji 



m, N B 

k,n o sc 



Moreover, define the following new random variable, Y k as: 

neq 

Thus, the achievable rate of the k th SS r k (t) is given by: 

D 

r k =-\og 2 (Y k ) (14) 

and equation (10) is reduced into 



A Y (s) = E{e'} 



(17) 



poo poo fO 

Jo Jo '"Jo 



00 S X-i Xr\ ...x 



1 -*2- A Af l 



X f Xl S X l)f Xl S X 2)-f Xl „ ( X M r )<VA 



x k,l 



^k,2 



k,M^ 



Define J t @X.X i+v ..X M and 

m k\~ l ~V^m 

A, I e 1 k,l e 



SXlh dx, 



(•00 jj; 

X @J jc/' 1 e v ' K +e" ri dx x (18) 

which is the innermost integral after dropping the constant 

1 / T{m kl )P k j ' , for simplicity. The integral in (18) can be 
rewritten in terms of the Meijer-G function [12] as: 

6> = f V« _1 Gjf [-£- rjG^i-sx, J 2 \-]d Xl (19) 



A 



A:,l 



Using the result provided in [13] for calculating 
integrals of hypergeometric type functions, one can find that 
the integral in (19) can be solved as: 
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n _ n m k,\ r U 



The integration on X 2 defined as 2 is given as: 



(20) 



^v*'* 2 ^ 



m k,2~ l ^~ X 2 , h,2 






"Ar.l 



dx, 



' m Jt,2" 1 "^^,2 



xG U - s fik,i J 3 X : 



3^2 



"Ar.l 



dx, 



Jo 



xg: 



X 



w fc,2 -1 ^l,0 r -^2 i-i 

2 Wl L"^ 'o J 

Pk,2 



S Pk, l^3 X 2 



1-ra 



fc,l 



rfx, 



(21) 

(22) 

(23) 

(24) 

(25) 



Similarly, using the result from [13], the integral in (21) is 
given by: 



#2 - Pk k \ l Pki 2 G d 



S Pk,lPk,2 3 



(26) 



Following the same argument used for calculating the second 
integral, the M-fold integral in (17), 6 M , is given by 



(27) 



1-ra, i ,...,l-m, ,, 
k,l k,M, 



inductio 


nas: 
f 

n 

K nGC k 


) 




xG 


,M k ,l 
l,M k 




A* 



Given # M , the closed-form formula of the MGF A F (s) 

k k 



is: 



A r » = - 







M,, 



(28) 



Substituting (27) in (28) we get the result in (16). 
The probability density function (pdf) of the random variable 

Y k is given in the following corollary. 
Corollary 1: The pdf of Y k is given by: 



(y-V-^X 



fy A (y) = - 



(y-i) 



Tift* 



m k,V m k,2'-' m k,M, 



T\ Y ^n) 

nsq 



(29) 
Proof: The pdf of Y k can be found as the inverse Laplace 

transform of A y (s) , that is 

(30) 



Using the formula for the inverse Laplace transform of the 
Meijer G-function from [14], we obtain the result in (29). 

Corollary 2: The cumulative distribution function (cdf), 

Fyfa (y) ' of Y k > is g iven b y ; 






Fy A (y): 



(y-D 



Tlftn 



m kx m k2 ,...,m kM 



n r K,j 



(31) 



Proof: By definition, the cdf F Y „ (y) is the probability that 
the random variable Y k is in the interval [ 1 , y] , that is 

F Y k K k ^ = \ y J Y A iz)dz (32) 



Therefore, 



'0,M. 









m k,V m k,2'-' m k,M, 



1 SloilrK^r) 



dz 



n r K,„)2^- ju !4 
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r(i-r) Y[p kr 



) T dr (33) 



neC 



where in the last step in (33) we used the equality 

r(l-r) = -zT(-r). 

U is an appropriate integration contour, and j = y—1 . 
Solving the integral in (??) will yield the solution in (31). 



calculated as the djfferenceJbetwepLthe minimum required 
data rate and4hjp attainable da.La rlt^Rfe^ltuftfted outage 
probability il Jpbialfied^ 4?&m (#4^ by substituting the 
bandwidth of n f& same SCs used in simulation per user and 
the Nakagami-m pdf parameters. Results are shown in Figure 
4 for the calculated outage probability against the simulated 
outage probability for different values of the minimum 
required data rates. The results in Figure 4 clearly show that 
the rate calculated outage probability matches that obtained 
from the simulation for different values of the user required 
data rate. Hence, the calculated rate outage probability is 
verified to provide a closed form expression of the rate outage 
probability in an OFDMA system. 

IV. Derivation of Multi-Hop Rate Outage 
Probability 



Using the CDF, F Y (y) = Pr(Y k < y) , in (31) 

k 

and (4), the closed form of the rate outage probability for a 
given Q is given by: 



pout = p 



** 



r 7 . IB 
/n k,min sc \ 



G M k ,l 



(2^" 



JB, 



(34) 



1) 



FK. 



"*,1''"*,2" 



n kM h 



n.C k 



n r Kj 



As an indicative example for the calculated outage 
probability, the outage probability is plotted in Figure 3 as a 
function of the minimum required data rate per user for 
different values of the channel physical bandwidth. The 
observed results evidently show that the outage probability 
decreases with the increase of the channel physical 
bandwidth for the same required data rate. Moreover, For 
larger channel band widths, for example when B is 10 MHz 
and 5 MHz, the difference in the outage probability values 
become tighter, as shown in Figure 3, with outage probability 

values of 2 x 10" 4 and 2 x 10~ 5 for a bandwidth of 5MHz 
and 10MHz respectively. 

We performed simulation experiments based on the 
system model described in Section 2 to compare the 
calculated rate outage probability in (34) against the rate 
outage probability evaluated from the simulation. The 
network consists of one BS and multiple SS. The BS 
allocates the available SCs per user based on the users' data 
rate requirement and their channel quality. This process is 
repeated each frame duration. For our simulation, we used a 
frame size of 10 ms. The matrix of SCs is generated every 
frame based on the Nakagami-m channel model and the SCs 
are allocated per user to meet the users minimum data rate 
requirements. The attainable data rate is calculated given the 
allocated SCs per user. The simulated outage probability is 



In this Section, using the results of section 3, we 
derive the rate outage probability of OFDMA multi-hop 
communication system under Nakagami-m channel models. 

For multi-hop scenario, a minimum requested 
transmission rate requirement for a given SS is satisfied when 
the minimum supported transmission rate on an intermediate 
link (bottle-neck) in the routing path between the BS and this 
particular SS is greater than the requested rate. Otherwise, a 
rate outage will occur. Therefore, the multi-hop rate outage 
probability is the probability that the supported transmission 
rate on the bottle-neck link is less than the minimum 
requested transmission rate by a particular SS. In order to 
mathematically express this probability, let us define the 
folio wings. For the intermediate link (hop) i , 

IG{ 1, 2,* * • 9 N h k } with N h k is the number of hops 

(distance) between the BS and kth SS, let 



IK 

"4 



(35) 



where Q is the subset of data subcarriers allocated to carry 
the data of the kth SS on the ith hop. Henceforth, the 
achievable rate per OFDMA symbol of the k th SS r[ on ith 
hop is: 

i=— iog 2 (y;) (36) 

The minimum supported transmission rate on the 
route from the BS to the SS k is defined as 
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r*"=min{rl,r?,:;r k h *} (37) 

From the above, we can figure out that the multi-hop rate 
outage probability of the kth SS, P™ m , is given by 

K U, - mh = PM n K<r k , min ) 08) 



or 



~)OUt-mh 



= p(Y k bn \q<2^ min,Bsc ) (39) 



Here, 



y h,k 



Q = IK 



Assuming that [Y k } i=1 ' are independent random variables, 
then 



-*out-mh 



= i-p r (r/ n i q > 2 rk ' min ' Bsc ) , 

N h,k , . 

= \-Y\p r \Y l k \q c >2 rk ' min sc \ 

i=\ 

N h k 

= i - ri f 1 - p r ( Y k ■ q < 2 rk ' min ' B$c )) 



(40) 



Note that P r y t 'iq<2 



.r,-IB, 



) 



is the rate outage 



probability at the ith hop, i.e. P™ ,l , thus (40) can be written 

as 



iy h,k 

p out-mh =l _Jj{l-P k ° Uti y 



(41) 



with 



l > M kJ +l 



T)OUt,l 

r k ~~ 



where M ki — 



r, ■ IB 
/r\ k,min sc 1 \ 



riA„ 

»4 



n k,V m k,2'-' m k,M 



k.i 



Yin 

«4 



m k J 



Q 



i.e. the number of data subcarriers 



allocated to carry the data of SS k on the ith hop. It should 
be noted that the assumption of the independence of 

\Xk h=i can resu lt f rom assigning each subcarrier once in 
the whole network, or from the multiuser diversity even if the 
same subcarrier has been assigned again in the network after 
some hops (to eliminate interference). 

In line with the single-hop simulation results 
provided in 4, we performed simulation experiment to 
compare the calculated multi-hop rate outage probability in 
(41). We used the same environment and simulation 
parameters as that used for the single hop. Our simulation is 
based on Type-II LTE-relay and the transparent 802. 16j relay 
mode. We chose these modes to allow centralized SCs 
allocation at the BS. Type-I LTE-relay and non-transparent 
802. 16j relay mode do not entail much difference from the 
point of view of verifying the correctness of the outage 
probability derivation. Note that we derived the closed form 
rate outage probability for a given user assuming the 
allocated sub-carriers' set is pre - allocated. The algorithm for 
allocating subcarriers to network users is out of the scope of 
this work. Algorithms for allocating subcarriers could be 
performed by the relay node in the non-transparent mode 
(distributed SCs allocation) or by the BS (centralized SCs 
allocation). For the sake of simplicity, we simulate here the 
transparent relay mode. Transparent relay mode restricts the 
number of hops in the network into 2 hops. The BS allocates 
SCs for the relay node and the SS given the Nakagami-m 
channel condition at each subcarrier at the relay and the 
access link. The attainable rate of the allocated SCs at each 
link is compared to the minimum required data rate to 
identify the bottleneck link. The simulated outage probability 
is calculated as the difference between the required data rate 
and the bottleneck attainable data rate. The calculated outage 
probability is calculated from (41) given the allocated SCs 
used in simulation and the Nakagami-m pdf parameters. 
Results are shown in Figure 5. Figure 5 shows similar results 
to those obtained for single-hop outage probability. The 
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calculated results almost match the simulated results for low 
and high required data rates. Hence, we can safely deduce 
that equation (41) provides a closed form expression of the 
rate outage probability in multi-hop OFDM A systems. 

V. OUTAGE PROBABILITY APPLICATIONS 

The closed form of the rate outage probability can be 
utilized in several aspects. For example, it can be integrated 
into a subcarrier allocation scheme, where the BS allocates the 
subcarriers and subcarriers' power level to the active SSs in the 
system to achieve a requested outage probability bound. Also, 
it can be used as a performance measure to evaluate the 
efficiency of different resource allocation algorithms, or to 
improve the performance of a certain resource allocation 
algorithm through adapting the resource allocation based on 
the feedback information of the channel status and the outage 
probability of users. As a final example, The SS rate outage 
probability can be used to calculate the network rate outage 
probability which can be used latter in efficient dynamic 
admission control algorithms. 

In this Section, as an example, we consider a 
centralized subcarrier allocation scheme that is aware of the 
outage probability of each user, called SCA-OUT, for single 
hop cellular network. In SCA-OUT, the BS allocates the 
subcarriers to the K active SSs in the system and then loads 
each SC with an optimal power level such that the overall DL 
achievable transmission rate is maximized. The constraint can 
be the exponential average transmission rate. This is assured to 
be bounded from above and from below by a given probability. 
Mathematically speaking, the problem is formulated as 
follows: 



**(') = ^Z(l-Vv,(n) (45) 



* n=0 * 



From (9), the achievable rate at the j th frame of the k th SS, 
r k (j) can be expressed as follows: 



r k (j) = ^log 2 (Y k (j)) 

= — iog 2 (i+r w U)) 



(46) 



(47) 



where Y Keff (j) = Y k (j)-l is the effective SNR over all k 

th SS's allocated subcarriers. Henceforth, / ke ff(j) is a 

random variable with the following probability density 
function: 



(yr l G^: 



Tin* 



'n kv m k2 ,...,m kM 



\/y > 0, (48) 

and accordingly, the CDF is given by: 



^l,M k+ l 



*wi (30= - 



y 



m k,V m k,2'-' m k,M k '° 



n r K,j 



,(49) 



SCA-OUT: "« ^ c ZI>g 2 (l + r M (0) (42) 



k=lriGC, 



s.t.: 



Pr(r kMn <R k (t)<p k r kmin \q)>s (43) 

Pk , n >o,q*0 

where l< P k < r kimvc / r kmin and 

R k (t) = (l-j)R k (t-l) + jr k (t), (44) 

is the exponential average received rate over a window size of 
T which is an integer multiple of the time frame T f . Using 

the fact that R k (— 1) = one can rewrite (44) as follows: 



To be able to solve the optimization problem in SCA- 
OUT, we need to provide a closed form expression of the 
probability given in (42). Therefore, we need to find the pdf of 

R k (t) , which is a weighted sum of independent random 
variables {r k (n)}"Z . Note that the independence of the 

random variables {r k (n) }^~ is a result of the assumption that 

the channel gains over all subcarriers between the k th SS and 
BS only change on a time frame basis. Next, we introduce 
some preliminary results that lead us to an approximate closed 
form formula of the constraint in (43). 

With the deployment of the AMC module at both the BS and 
SS to keep the bit error rate (BER) at a prescribed value, the 
achievable rate is a discrete random variable, see Table 1. 

Henceforth, the probability mass function (pmf) of r k {jl) at 
the n th time frame denoted to *¥•, and is given by: 
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x P i @P r (r k (n) = v.) (50) 

= P r (SNR t < Ykeff (n) < SNR M ) (51) 
C SNR i+i 

fy \c (z)dz 
isNR ( J r k ,effU) l( i 

where {SNR ( }. =1 are the thresholds of the 7 non-overlapping 
regions shown in Table 1 with y 8 = oo , y Q = , and 
V. E { 0.5, 1 , 1 .5, 2, 3, 4, 4.5 } is the number of information 

bits on the M -ary QAM symbol. Providing a closed formula 
for the probability in (42) requires the probability density 

function (pdf) of R k (t) . Assuming that first and second order 
moments of R k (t) are finite, for a large t , and according to 

the Central Limit Theorem (CLT), the pdf of R k (t) can be 

approximated as a Normal distributed random variable as 
follows: 

*11->N(0,1), (52) 

a R k (t) 

where N (0, 1) denotes a normal distribution with zero mean 

and variance equal to 1 , with JU R (t) being the mean of 

R k (t) and Or t t \ is the variance of R k (t) . Both jU R (/) and 



'Rk(t) 



can be found as functions of Li , , and <T ,. which 



are the mean and the variance of T k {n) , respectively: 
1 t 1 

1 n=0 l 



(53) 



i=\ 



T n =0 * 



and 



°V> " 



1) y(l_l)2(r- n) 2 



n=0 



where 



<(«)=^K 2 W}-^(,) 



(54) 



^ 7 



i=l V ?=1 

Therefore, the approximate pdf of the random variable 
R k (t) is given as: 



/*(o (*) = 



v^Vrm 



-exp 






(55) 



Next we rewrite an approximate closed form expression for 
equation (42) as follows: 

Pr(r kMn <R k (t)<p k r kMn \Q) 

?k k,min 



= [ K """" f RtM (z)dz 

k,min 



f ( 
erf 

\ v 



P r k,min ^ R k (0 



4^cj r 



-erf 



(56) 



f _ \\ 

r k,min f^R k (t) 



2 



^RAt) 



JJ 



error 



2 f© 2 

where erf(&) = — i=\ exp(-x )dx is the 

v^ ° 

function. Then, the problem SCA-OUT can be rewritten as 
follows: 



SCA-OUT:^ J B s X^og 2 (l + r k , n (t)) \ (57) 



k=ln&Ci 



s.t.: 1 

2 



X e/ / 



Ph. 



' ^R k 0) 



V V 
f 



V2~<x s , 



(0 



•erf 



*& a RM 



>£ 



J J 



Pk >o,q*0 



(58) 



This problem can be either solved analytically or by a heuristic 
algorithm, which is beyond the scope of the paper. 

VI. Conclusions 

We introduced an analytical model to calculate a 
tractable closed form rate outage probability expression of 
a SS in a single and multi-hop cellular network. The 
closed form expression of the rate outage probability is 
obtained by analyzing the desired data rate of each user 
against the link attainable transmission rate. The problem 
is formulated as a function of SNR per OFDMA subcarrier 
and the Nakagami-m channel model. Simulation results 
verified the match of the calculated outage probability 
with the simulated one in an OFDMA system. Moreover, 
we formulated an optimization problem for single -hop 
scenarios in which subcarriers are assigned to admitted 
users in the network such that the total transmission rate is 
maximized while catering for fairness for all users. 
Fairness was achieved by guaranteeing a minimum 
average rate outage probability over a time window of 
length equal to integer multiple of the time frame for each 
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admitted user. 
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Table 1: Modulation and Coding Schemes for IEEE 802.16. 



Modulation 
(coding) 


Info 
bits/symbol 


Required SNR 


BPSK(l/2) 


0.5 


6.4 


QPSK(l/2) 


1 


9.4 


QPSK(3/4) 


1.5 


11.2 


QAM(l/2) 


2 


16.4 


QAM(3/4) 


3 


18.2 


QAM(2/3) 


4 


22.7 


QAM(3/4) 


4.5 


24.4 
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Figure 4: Single hop Outage Probability against 
the Required Data Rate, calculated and simulated 



Figure 1 : Single hop cellular radio network 
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Figure 2: Multi-hop cellular radio network abstraction 



Figure 5: Multi-hop Outage Probability against 
the Required Data Rate, calculated and simulated 
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Figure 3: Single hop Outage Probability against the Required 
Data Rate for different channel bandwidth 
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